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ABSTRACT 

The  model  of  gravitational  surface  waves  on  beaches  of  small  slope 
formulated  in  Part  I  and  its  mathematical  theory  are  used  to  show  how  an 
incident-wave  amplitude  can  be  defined  so  that  a  bound  on  it  guar tees 
solutions  which  respect  the  assumptions  of  the  model  everywhere  and  forever. 
The  structure  of  those  solutions  “far"  from  shore  is  then  compared  with  that 
predicted  "near"  shore  by  the  classical,  linear  theory’  to  remove  the 
indeterminacies  of  both  theories:  Shore  reflection  is  determined  for  the 
classical  theory,  and  it  is  shown  how  the  critical  length  scale  and  amplitude 


•  •  •  •  • 

•  .•  V  J 


r. 

t-  •->.  : 


of  the  beach  theory  are  related  to  the  familiar  wavelength  and  amplitude  in 
deep  water.  These  results  indicate  that  the  beach  theory ^ captures  and 


elucidates  the  basic  singularity  structure  underlying  the  shore  behavior  of 
gravitational  surface  waves.  ^  J-'y.'t-..  ,  . 
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SIGNIFICANCE  AND  EXPLANATION 


The  shore  processes  of  oceanography  are  of  decisive  Importance  to  the 
containment  of  beach  erosion  by  the  Corps  of  Engineers*  They  are  also  of 
interest  to  Geology  in  regard  to  continental  accretion  by  sedimentation,  and 
to  Biology,  in  regard  to  the  ecology  of  coastal  wetlands  and  to  the  origins  of 
life  at  the  land-sea  interface. 

A  question  central  to  all  those  matters  is  how  water  waves  behave  on 
natural  beaches  close  to  shore,  and  progress  in  coastal  oceanography  has  been 
long  held  up  by  a  lack  of  theoretical  understanding  of  that  wave  mechanism. 
Waves  on  beaches,  of  course,  have  attracted  mathematical  Interest  for  a  long 
time,  and  it  has  been  recognized  gradually  that  the  shore  singularities  of  the 
classical  theories  are  an  incorrect  description  of  the  shore  mechanism. 

The  present  work  aims  to  make  a  convincing  contribution  to  the 
identification  and  elucidation  of  the  real  singularity  which  embodies  the 
underlying  structure  of  water  motion  close  to  shore.  It  is  an  essentially 
nonlinear  mechanism,  but  one  very  different  from  the  nonlinearity  of  waves  of 
large  amplitude  far  from  shores.  It  dominates  waves  of  small  amplitude  as 
well,  and  this  report  concentrates  on  unbroken  waves  to  complement  earlier 
explorations  of  breaker  collapse  and  conversion  to  run-up  to  a  clear  and 
balanced  picture. 

For  unbroken  waves,  the  new  beach  theory  developed  in  the  Reports  TSR  #'s 
2871  and  2872  is  especially  successful  because  it  will  here  be  shown  to 
overlap  with  the  classical  theory  of  small-amplitude  waves.  That  produces  a 
complete  mathematical  description  of  the  inviscid  structure  underlying 
unbroken  waves  over  natural  beaches,  which  reaches  from  deep  water  to  the 
moving  shoreline.  It  encompasses  standing  waves  as  well  as  waves  developing 
almost  arbitrarily  in  time,  and  shows  plainly  how  the  nonlinear  shore 
singularity,  which  manifests  itself  fully  only  close  to  the  moving  shoreline, 
nevertheless  determines  wave  structure  far  into  the  deeper  water. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


OM  THE  SHORE  SINGULARITY  OF  WATER-WAVES  THEORY. 
PART  II.  SMALL  WAVES  DON’T  BREAK  ON  GENTLE  BEACHES 


R.  E.  Meyer 

I.  Introduction 

A  model  for  gravitational  surface  waves  close  to  shore  over  beaches  of  small  slope  was 
formulated  In  Part  l'  (hereafter  referred  to  as  [F.  1  with  note  of  section  or  equation 

number).  Ihat  formulation  allows  for  the  study  of  waves  which  "do  not  brea)c",  that  Is,  of 


solution  of  the  nonlinear  beach  equations 


K  *  I;  ""■>  ■  “ 


(1) 


(2) 


which  satisfy  the  premises  of  the  model  [F.  IV]  In  a  strict  sense.  The  notation  and  Its 

physical  connotations  are  explained  at  length  In  (F.  IV,  V] .  The  mathematical  theory  of 

2 

the  model  Is  developed  in  a  related  report  (hereafter  referred  to  as  [M.  ]).  The 

account  to  follow  now  concentrates  on  the  oceanographical  Interpretation  and  application  of 
the  mathematical  theory. 

After  a  general  characterization  of  the  main  class  of  model-preserving  solutions  in 
Sections  II,  III,  the  account  turns  to  the  basic  of  removing  a  serious  shortcoming  of  the 
theory:  The  general  model  has  physical  validity  only  In  a  close  neighborhood  of  the  moving 
shoreline,  the  extent  of  which  It  cannot  determine.  That  is  a  fruitful  restriction,  on  the 
one  hand,  because  it  directs  the  analysis  toward  the  Key  to  wave  structure  on  beaches  of 

Q 

small  slope.  It  also  recognizes  that  the  model  does  not  describe  the  observed  mechanisms 
of  wave-breaking  during  shoaling.  It  is  a  great  handicap,  on  the  other  hand,  because  It 
leaves  unclear  how  the  input  to  the  model  could  be  determined  from  what  is  known  or 
observable  of  waves.  A  major  objective  now  will  be  to  show  that  this  difficulty  is  not 
present  for  model-preserving  solutions,  on  the  contrary,  the  model  describes  for  them  the 
structure  of  most  of  the  shoaling  process. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041  and  partly 
supported  by  the  National  Science  Foundation  under  Grant  MCS-821S064  and  a  Senior 
Fellowship  of  the  British  Science  and  Engineering  Research  Council;  the  author  is  also 
indebted  to  the  Imperial  Ck>llege,  London,  for  hospitality  during  the  drafting  of  the 
report . 


To  explain  this.  It  Is  advantageous  to  restrict  attention  In  Sections  IV,  V  to  the 
special  case  of  standing  waves,  for  which  the  discussion  can  be  concrete  and  explicit. 
There,  the  extension  of  the  model-preserving,  standing-wave  solution  beyond  the  model's  a- 
prlorl  domain  of  validity  Is  compared  with  the  extension  of  the  solution^  of  the  classical 
linear  theory  of  small-euaplitude  waves  on  beaches  to  the  smallest  shore  distances  at  wl  .ch 
that  theory  can  describe  water  waves.  The  two  solutions  turn  out  to  have  the  same  fl' 

(and  close)  approximation  over  an  Interval  of  distances  If,  and  only  If,  a  certain 
condition  (Equation  (17)  below)  Is  satisfied. 

This  condition  removes  the  basic  Indetermlnacles  of  ]cx>th  theories  by  determlnlm  n 
the  one  hatid,  the  correct  linear  combination  of  the  fundeunental  solutions  of  the  lines 
theory,  and  on  the  other  hand,  the  unknown  length  scale  L  of  the  beach  theory  IP.  IV]  . id 
the  beach  ^unplitude  (Section  II)  in  terms  of  the  deep-water  wavelength  and  amplitude.  A 
full  description  of  the  model-preserving  standing  waves  from  the  deep  sea  to  the  moving 
shoreline  Is  thereby  established.  It  Is  used  In  Section  V  to  derive  also  some  simple  and 
illuminating,  if  rough,  quantitative  predictions  about  the  amplification  of  such  waves  by 
shoaling.  Altogether,  standing  waves  are  simple  and  specific  enough  to  display  a  good  deal 
of  physical  flesh  (Sections  IV,  V)  on  the  powerful,  but  somewhat  abstract,  mathematical 
bones • 

Section  VI  returns  to  the  discussion  of  waves  of  much  more  general  time-dependence. 

For  model-preserving  solutions,  the  "Friedrichs  Connection"  between  beach  theory  and 
classical,  linear  theory  is  there  shown  to  cover  the  case  of  quite  general  time-dependence 
to  the  extent  that  the  waves  have  a  Fourier  transform  with  respect  to  time  in  deeper 
water.  Between  them,  the  two  models  offer  a  complete  shoaling  description  also  for  such 
general  waves  developing  in  time.  The  condition  for  model-preservation  (Section  VI)  is 
then  mainly  a  restriction  on  the  high-frequency  part  of  the  amplitude  spectrum. 


I 1 .  Ampl itude 

If  a  solution  of  the  beach  equations  (1),  (2)  Is  to  preserve  their  physical  validity 
as  an  approximate  description  of  water  waves  on  beaches,  it  is  necessary  that  the 
associated,  apparent  solution  be  Invertible,  as  will  be  seen  clearly  in  the  next  Section. 
For  this,  in  turn,  the  Inversion  Criterion  (M.  S,  F.  X)  states  a  necessary  and  sufficient 
condition  for  the  class  of  incident  waves  specified  by  3-compatlble  data: 

11m  8t/3X 

a*0 

must  have  a  positive  lower  bound.  It  is  a  rigorous,  but  not  a  practical  condition.  A 
useful  one  should  be  in  terms  of  a  relatively  simple  property  of  the  incident  wave.  Of 
course,  the  Riemann  representation  established  in  the  existence  proof  [M.  3]  offers  a  means 
of  representing  lim  at/3X  as  a  functional  of  the  data  a(a)  specifying  the  Incident 
wave,  but  the  resulting,  lengthy  formulae  are  too  complicated  for  practical  use.  Those 
data,  moreover,  specify  the  incident  wave  at  an  unknown,  and  possibly  quite  small,  distance 
from  shore.  They  need  to  be  related  to  oceanographically  meaningful  properties  at 
distances  where  the  Incident  wave  can  be  reasonably  observed  or  generated. 

In  selecting  a  practical  amplitude  concept  for  the  characterisation  of  model¬ 
preserving  solutions,  it  is  important  to  understand  the  choices  and  limitations.  An 
Important  limitation  arises  from  the  fact  that,  for  a  nonlinear,  singular  wave  equation 
such  as  (1),  (2),  the  eunplitude  concept  depends  on  the  solution  class.  It  is  different, 
for  instance,  for  the  highly  singular  solutions^' ^  related  to  "bore  collapse”  and  for  the 
less  singular  solutions^*'*  for  (-1 ) -compatible  data:  A  correct  2unplitude  concept  is 
essentially  hinged  to  the  fine  print  in  its  definition.  The  mathematical  theory  (M.  5] 
shows  plainly  that  1-compatibility  is  insufficient  for  global  inversion  of  apparent 
solutions,  and  no  general  criterion  for  it  has  yet  been  found  for  2-comparible  data.  At 
present,  a  reliable  amplitude  concept  for  model-preserving  solutions  must  therefore  be 
restricted  to  the  class  of  3-compatible  incidence  data.  This  is  quite  adequate  for 
oceanographical  purposes  because  it  requires  [M.  3]  much  less  smoothness  than  waves  are 
expected  to  preserve  during  shoaling  short  of  breaking. 


A  related  restriction  is  that  the  data  auat  be  consistent  with  local  inversion  on  the 


incidence  boundary.  For  that 


(3) 


is  necessary  and  sufficient.  Sufficient,  because  the  Jacobian  J  ••  3(x,t/3(a,$)  of  the 


characteristic  transfomatlon  (F.  V]  is 


J  -  j  a((3t/3X)^  -  (3t/3c)^J 


(4) 


by  [F.  (4),  (6),  (12)],  and  Ijy  [F.  (8),  (12),  (13),  (20)1,  3t/3<J  -  0,  3t/3X  -  1/2 

initially;  by  (F.  (8),  (12)],  (3)  leaves  3t/3o  +  3t/3X  >  0,  and  by  [F.  (8),  (9),  (12)], 
3t/3a  -  3t/3X  then  remains  negative  along  the  incidence  boundary  0  -  X  -  const  ■*  2a.  The 


1  -3/2  * 

inequality  (3)  is  also  necessary  because  3t/3o  3t/3X  «  —  +  2o  a  would  otherwise 


have  a  root,  and  by  [F.  (8),  (9),  (12)],  3t/30  —  3t/3X  would  still  be  negative  at  the 

first  such  root  and  it  would  be  a  root  also  of  the  Jacobian. 

Among  the  choices  is  whether  amplitude  should  refer  to  velocity  or  surface  elevation 


or,  as  has  come  to  be  preferred  often  In  oceanography,  to  an  intrinsically  ncndisansional 
quantity  such  as  surface  slope  3r\/3x.  The  last  comends  Itself  by  the  Intimate  connection 
betvraen  3n/3x  and  fluid  acceleration  documented  by  the  momentum  balance  (2)  and  will  be 
adopted  here.  That  prompts,  in  turn,  that  amplitude  be  formulated  as  a  property  of  the 
acceleration  measure  a  characteristic  of  the  incident  wave  [F.  V]  and  in  particular,  of 
the  function  a(2o-2a)  in  [F.  (22)]  which  specifies  the  Incident  wave  in  the  mathematical 

beach  theory  [H.  IX] .  The  simplest  choice  is  then  the  following. 

*  —  3/2  — 

Definition.  For  3-compatible  incidence  data  a(a)  >  -(ato)  /4  on  [0,a] ,  beach 


amplitude  means 


'  »  •  ^  . 


III.  Small  waves  never  break 

Inversion  Ibeorem.  An  apparent  solution  is  globally  invertible  if  the  beach  amplitude 
5  is  not  too  large. 

This  is  a  direct  corollary  [N.  6]  of  the  Inversion  Criterion  because  the  canonical 

equations  [F.  (9),  (10)1  are  linear  and  homogeneous,  so  that  the  characteristic 

accelerations  a  and  b  scale  in  proportion  to  S  throughout  the  apparent  domain  [F. 

(23)1,  and  because  [F.  (13)1  and  the  Regularity  Corollary  R. 3a  [M.  4]  show 
1  “3/2 

3t/31  -  —  ”  d  (a-b)  to  have  a  limit  as  o  ♦  0,  which  must  also  scale  with  6. 

In  the  first  place,  the  Inversion  Theorem  documents  only  the  opportunity  for  a  direct 
examination  of  the  apparent  solutions  in  terms  of  the  beach  equations  (1),  (2):  the 
functions  x(a,6)  and  t((>,0)  obtained  by  quadratures  from  Y(a,X)  and  Z(a,X)  via  [F. 
(13),  (9),  (10),  (12)  and  (6)1  then  furnish  a  regular,  pareunetrlc  representation  of 


velocity 

u( x,t)  =  ■“  ( X+o)  -  t 

and  surface  elevation 


2 

n(x,t)  “  h  +■  x  ”  o  /16  +  X  . 

Beyond  this,  however,  the  inversion  condition  is  a  mathematical  code-word  for  a 
physical  bound.  By  (F.  (11),  (12)]  and  (4),  (2),  the  surface  slope  is 

3n/3x  »  1  -0(3t/3X)/(4J) 


(5) 


”  -3u/3t  -  u3u/3x  . 

To  find  a  root  of  the  Jacobian  J  therefore  means  that  the  beach  equations  (1),  (2)  have 
already  broken  down  as  a  physical  model  [F.  IV) :  invertibllity  is  a  necessary  condition 
for  preservation  of  the  model.  Small  beach  amplitude  (  assures  that  this  necessary 
condition  remains  satisfied  permanently. 

It  is  not  a  sufficient  condition,  however.  The  real  issue  must  be  whether  fluid 
accelerations,  surface  slopes  or  surface  curvatures  occur  which  are  quantitatively  large 
enough  to  shake  our  confidence  in  the  scaling  (F.  IV)  that  leads  to  the  beach  equations  as 
an  approximate  description  of  water-wave  properties.  Ihysical  effects  ignored  in  the 
derivation  of  this  model  must  then  be  considered,  and  they  must  determine  the  interpreta- 


-5- 


tion  of  the  phrase  "quantitatively  large  enough."  Any  sufficient  condition  for  confidence 
in  the  model  must  therefore  involve  quantitative  bounds  on  the  surface  slope. 

Since  the  Invertibility  Corollary  [F.  X]  shows  the  threat  to  the  model  to  be  greatest 

near  the  shoreline,  it  is  appropriate  to  look  for  surface-slope  bounds  there,  and 

-3/2  1  -3/2 

lim  3n/3x  “  lim  [a  ® 

0*0  o+O 

=  -lim  (3u/3X)/(3t/3X)  , 

0*0 

by  (5),  (4),  [F.  (12),  (13)1  emd  the  Regularity  Corollary  R3  [M.  IV) .  The  righthand  limit 
exists,  by  the  Inversion  Criterion,  and  since  Y  scales  in  proportion  to  5,  the  surface 
slope  has  a  permanent  t)Ound  even  at  the  shoreline.  Hence,  if  physical  considerations 
indicate  a  tolerance  limit  for  surface  slope  magnitude,  then  a  positive  bound  on  the 
amplitude  6  exists  which  will  guarantee  that  the  tolerance  limit  is  respected. 

Again,  the  theory  [M]  permits  that  amplitude  bound  to  be  evaluated  from  the  Rlemann 
representation  of  the  apparent  solution,  but  a  very  general  prediction  of  the  bound  is  too 
complicated  to  promise  much  profit.  It  must  also  be  borne  in  mind  that  the  model  of  a 
beach  of  uniform  slope  is  designed  to  serve  conceptual  clarification  more  than  direct 
description  of  real  circumstances.  Moreover,  amplitude  bounds  on  the  incident  wave  in  the 
close  vicinity  of  the  shore  covered  by  the  theory,  up  to  this  point,  are  not  of  direct 
value  to  oceanography.  For  these  reasons,  attention  will  be  directed  now  to  the  prediction 
of  bounds  in  more  specific  circumstances  of  practical  interest,  and  to  the  extension  of 
such  predictions  to  bounds  on  the  incident  wave  in  deeper  water,  which  are  of  oceanographic 
value.  To  this  end,  it  is  of  great  help  to  consider  first  a  very  special  exeunple  admitting 
a  concrete  elucidation  of  such  issues,  before  the  discussion  returns  in  Section  VI  to  a 
much  wider  subset  of  the  very  general  class  of  wave  motions  encompassed  by  the  mathematical 


theory. 


IV.  standing  Waves 


Carrier  and  Greenspan's  Iheom'^.  If  the  beach  equations  (1),  (2)  have  an  apparent 


solution  with  velocity  simple-harmonic  in  characteristic  time,  its  velocity  u  and  surface 


elevation  n  must  be 


u  =  o“’3^/3o,  n  -  jctoVie  “  —  34i/3\  -  -  u^ 

4  2 


respectively#  with  X«2t  +  2u-<J/  a*4h 


♦  -  85 


and  constants  d  >  0#  6  >  0  and  o)e  For 

<  1  , 

moreover,  this  apparent  solution  does  represent  a  strong  solution  u(x,t),  h(x,t)  of  the 

3  3 

beach  equations  (1),  (2),  which  has  a  bound  |ai  d|/(1-{<ii  6|)  on  the  magnitude  of  surface 
slope  and  fluid  acceleration. 

Such  a  standing  wave,  of  course,  can  be  realized  only  asymptotically  as  t  +  *,  and 
the  present  theory  (Section  II,  III)  adds  tiic  coBunent  that  it  can  develop  in  time  frost  the 
undisturbed  state,  if  the  incident  wave  amplitude  remains  small  enough  during  the 
development,  and  that  it  must  then  develop,  if  the  incident  wave  tends  ultimately  to  the 


'i 


incident  part  of  the  standing  wave. 

The  proof  implicit  in  the  original  account^^  will  help  both  to  illustrate  the 
application  of  the  Inversion  Criterion  and  to  prepare  the  discussion  to  follow.  In  terms 
of  the  independent  variables  a,  X  the  characteristic  form  [F.  (6)]  of  the  beach  equations 


3x  3t  V2 
3o  “  “  3o  “  "  31' 


V2  3t  3t 
Ja  *  Tx  ’ 


and  if  functions  9<^/3o  and  d^/9X  are  defined  in  terms  of  u  and  x  as  in  the  theorem, 
the  last  two  equations  are  equivalent  to 

2 

9  9^  9  9<^  9  9<}>  9 

3a  3X  '  3X  3o'  lo  lo’  “  ‘ 

If  u  is  to  be  simple-harmonic  in  X  with  frequency  u)/2,  say,  it  follows  from  the  last 


equation  that 


ii(o,X)  =  Zg('VO/2)e 


-iU)X/2 


-7- 


In  terms  of  Bessel  functions  or  order  n, 


Z  (^)  -  8«  J  (|2.)  +  84  Y 

n  «  n  2  2  n  2 


and  from  the  definition  of  3^/3<j,  3^/31,  that 


“  »  <<00.  -i<iiX/2 


ioi  _  ,(iK>.  -i<i)X/2 
8-  ^O'— 


2t  =  X  +  0-  2u-X  +  o+  —  Z,  . 

0  12 

If  ^2  ^  this  apparent  solution  is  not  admissible  (F.  X]  because  clock-time  t  does 
not  then  exist  at  the  shore  a  =  0. 

With  $2  °  however,  the  solution  la  admissible  (F.  X],  and 

1 


3t 

3X 


216(11  .aw,  -i(i)X/2 

- - J  ( — )e 

2  a  12 


Since  is  then  analytic  in  a  and  X,  the  incidence  data  eure  n-compatible  for  all 

n  >  0,  and  the  Inversion  Criterion  (Section  II)  becomes  here 

lim  3t/3X  “  4  7 

a*0  ^  ^ 

for  some  6^  >  0,  which  is  satisfied  if,  and  only  if,  For  this  simple, 

apparent  solution,  it  can  also  be  checked  readily  that  the  Jacobian  (4)  has  roots  at 

positive  values  of  o  if,  and  only  if,  |w^6|  >  1. 

By  (4)  and  (5),  moreover,  the  surface  slope  is 

3n/3x  -  1  +  -J  (3t/3X)/((3t/3a)^  -  (3t/3X)^] 

and  from  the  standard  properties  of  Bessel  functions  it  is  seen  to  have  magnitude 

|3n/3x|  <  sup  lim  |3ri/3xj  =  |(i>^6|/(1  -  lii)^5|)  , 

X  0*0 

so  that  the  bound  on  6  corresponding  to  any  desired  bound  on  l3r|/3x|  is  readily 
determined.  (It  will  emerge  in  the  next  Section  that  there  is  no  loss  of  generality  in 

normalizing  the  non-dimensional  frequency  ui  to  1  for  waves  of  fixed  frequency.) 

1  2 

Naturally,  this  standing-wave  solution  has  many  interesting  properties  ,  of  which  the 
most  important,  perhaps,  is  the  asymptotic  decay  with  increasing  distance  from  shore.  From 


(8) 


the  asymptotics  of  the  Bessel  functions  as  4|(»>|h 


=  |(olo 


the  explicit  formulae 


just  quoted  are  seen  to  imply,  i*.  turn 


“3/2  ,Ei>o  lit  -x(uA/2 

-86(— )  ^  0  cos{- - )  e 

ir  2  4 


n 


>2i5( 


UjV2 


-  V2  , 

0  cos(— “ 


IT  -iujX/2 

v)  e 


uiX/2  ~  (i)t  -  UO/2  + 


|x|  =  O^/ie  +  0(0  (1IO/2  =  2<ii|x|^^+  0(|x|  ,  (9) 

u  ~  |x|  cos(2;o|x  -  3it/4)  e  . 

n  ~  |xl  cos(2u|xl^-  it/4)  e  •  ('0) 

The  motion  far  from  shore  (in  units  of  the  still  unknown,  critical  length  scale  L)  there 
fore  becomes  simple-harmonic  also  in  clock-time  and  the  amplitude  of  the  surface  elevation 
comes  to  conform  to  the  Green's  law  typical  of  linear  longwave  (or  tidal)  theory. 


V.  Friedrichs  Connection 


Ihe  beach  equations  (1),  (2)  etre  in  terms  of  the  nondimensional  variables 

•  *  •  > 

«  X  /X,  tg  -  Ut  /X  , 

h„  -  h*/(eL)  -  0^/16,  n„  -  h„  +  x„  -  n*/(eL)  J 

B  B  B  B 

•  2 

=  (iS^t  /tg  »  Uil^/V,  U  “  egL 


where  stars  denote  the  dimensional  quantities,  and  the  critical  length  scale  L  still 
remains  undetermined.  The  seaward  asymptotics  of  the  beach  theory  describes  the  solution 
character  towards  somewhat  deeper  water  at  somewhat  larger  distances  from  shore  in  the 


limit  m  o/2  +  «.  when  the  local  2Lmplitude  decreases  with  distance  from  shore,  as  for 
B 


standing  waves  (Section  IV),  e«ga,  then  ^  the  limit  is  also  described  by 


B 


(ij  a 

B 


“.L  /1V2 

~  2  -  — I  ♦  “  • 

U  '  ' 
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In  dimensional  notation,  the  first  approximation  to  the  surface  elevation  (6)  of  the 
Carrier-Greenspan^^  standing  wave  is  then 


-i«eL  ( 


2L  ,V2 


IT' 


eX 


1V4 


2irx 


’(2 


2iix 

eX_ 


iV2 


’ll 

-  e 


-iu^t 
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by  (9 )  (11),  where 


X^  =  2irg/<fl, 


is  the  dimensional  "deep-sea"  wavelength. 


The  solution  of  the  classical,  linear  theory  of  irrotational,  gravity  waves  over  a 

3 


beach  of  uniform,  small  slope  e  has  been  given  by  Friedrichs  .  His  notation  uses  the 
variables 


(Xp,Zp)  =  (2ii/X^)(x  ,z  ), 


tp  =  -m.t 


n  =■  2xn  /X_ 
F  « 


and  like  the  beach  theory,  his  "second  asymptotic"  description  explores  the  limit  E  ♦  0 

•  • 

with  Ex  fixed,  but  at  a  possibly  different  order  of  ex  ,  namely,  on  the  scale 


X The  regular  standing  wave  of  the  classical,  linear  theory  is  then  described  by 


V  ->  ...  1  .n.  L  V  A  LU.M  gl-Jl.'  >.1 J,  J  J  J I .  ■  J 1.  U .  IJ  L .  j  i  y  ... .., 


It, 


-  21«j^a^(A)  coa(^^  -  j)  e  ^  (14) 

In  terms  of  certain  functions^  of  the  local  wavelength  1  and  of  a  normalization  factor 


AV-  A. 
•'»J 


Wie  "deep-sea"  asymptotics <  as  ex^  ♦  *,  of  this  description  Is  given*’  by 


(2it£ 


h. 


ex^  t  r  /8  and  therefore,  by 


(15) 


-  -21«jj(2tc)‘^2  8in(xp  -  Ij'  -  f )  e  ^  “■  "  > 

so  that  the  "deep-sea*  wavelength  Is  just  and  the  "deep-sea"  amplitude, 

A  diversion  seems  necessary  here  to  clarify  the  role  of  the  "deep-sea  limit*  in  the 
theory  of  waves  on  beaches.  It  is  not  thought  of  as  a  realized  wave  state,  but  as  the  most 
natural  and  convenient  reference  state  Independent  of  the  circumstances  of  a  particular 
experiment  or  particular  offshore  topography.  To  the  extent  that  the  classical,  linear 
theory,  and  camnon  approximations^^  for  naturally  gentle  sealsed  topographies  can  serve,  a 
simple  algorithm  —  which  Is  the  more  accurate^,  the  smaller  the  beach  slope  ~  relates 

this  reference  state  to  the  local  wave  state  at  the  position  on  a  beach  where  the 

* 

undisturbed  depth  is  ex  .  This  Is  the  sense  In  which  the  deep-sea  values  of  wavelength 
and  amplitude  for  a  beach  of  uniform  slope  were  used  as  reference  by  Friedrichs^  and  will 
continue  to  be  used  here. 

Friedrichs*  account  makes  clear  why  the  scaling  underlying  the  classical,  linear 
theory  precludes  its  extension  to  the  shore.  The  description  (14),  accordingly,  applies 
only  for  cx^,  >  const  >  0,  but  there  Is  no  mathematical  obstacle  to  the  choice  of  a  small 
number  for  this  constant  and  to  the  exploration  of  (14)  for  values  of  cx^  near  such  a 
number.  The  first  approximation  to  (14)  for  values  of  sx^  small  in  this  sense  is^ 

it 


2  Va  X  Vj  ic_ 

R  o,  X  ct  e>  '»  f,,  F,  ^  X-.  F 

~  2i5„( - J  cos(2(— )  -  tJ  e 

F  R'x  e  A' 

r 


,v. 


-iu>,t 


(16) 
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«.•  *✓  V  *J! 


Nr 


m 


V' 


V.-.V.  .o 


<17) 


*  *  * 

which  represents  the  same  function  ti  (x  ,t  )  “  as  (13)  If  and  only  If, 

-  «,/(e«)  • 

Together  with  the  amplitude  bound  of  Carrier's  Theorem  (Section  IV),  (17)  therefore  gives 
the  necessary  and  sufficient  condition  for  a  coherent  and  rigorously  founded  approxlsiatlon 
to  the  description  of  unbroken,  gravitational,  standing  waves  without  vortlclty  over  a 
plane  beach  from  the  deep  sea  right  to  the  moving  shoreline. 

It  Is  disappointing,  but  also  Illuminating,  that  (17)  establishes  only  a  relation 

between  the  scale  L  and  the  shore  amplitude  S,  given  the  reference  wave  farther  from 

shore.  The  beach  model  Is  more  subtle  than  could  have  been  suspected  at  the  beginning 
[F,  IV],  It  Involves  an  organic  connection  between  length  scale  and  euiplltude. 

This  does  not,  however,  d2Uiiage  the  determinacy  of  the  model  In  regard  to  observables. 
For  Instance,  If  amplitude  Is  thought  of  In  terms  of  energy,  then  It  may  be  measured,  far 
from  shore,  by  the  maximum  of  the  surface  elevation  |n  |  so  that  the  deep-sea  amplitude 

near  shore.  Carrier  and  Greenspan's  Theorem  (Section  IV)  shows  the 

correct  measure  to  be  the  maximum  of  |n  -Ki  /(2g)|,  %«hlch  the  theorem  shows  to  be 
EUOgfi,  l,e., 

Vo 

(V(2e)]'2 

-  Vo 

times  the  deep-sea  amplitude,  by  (17),  The  amplification  factor  e  for  beach  slopes 
has  been  Identified  also  In  more  heuristic  descriptions^^.  Similarly,  the  Identification 
of  (13)  with  (16)  by  means  of  (17)  demonstrates  how  knowledge  of  frequency  and  reference 
amplitude  suffices  to  determine  the  detailed  shore  behavior  of  unbroken  standing  waves 

completely  via  Carrier  and  Greenspan's  Theorem  (Section  IV), 

An  even  more  rational  definition  of  amplitude  Is  In  terms  of  the  non-dlmenslonal 

*  * 

surface  slope  3n  /3x  ,  which  Is  also  of  prime  dynamical  relevemce,  by  (2),  In  the  deep 
sea,  meuc|3n  /3x  [  Is  just  1 5^1  ,  and  at  the  shore.  It  Is  e|ii)^5|/(1  -  |ii)g5|),  by  (9); 
and  by  ( 17 ) , 


«  (it/2)^2e“5/2  j 

B  “ 

-3/2 

The  eunpllf icatlon  of  surface  slope  Is  therefore  e  /1(2/ti) 


,V2.  -5/2 


5  ] ,  and  the 


inversion  condition  l“g^l  ^  ^  necessary  and  sufficient  for  unbroken  standing  waves 


o 

translates  Into  the  deep-sea  condition 

This  expresses  the  true  sense  in  which  the  amplitude  must  be  small  for  waves  that  "never 
break"  [F.  II]:  the  amplitude  far  from  shore  must  be  small,  not  the  shore  amplitude 
The  precise  condition  just  given  is  unrealistically  stringent,  however,  because  viscosity 
and  surface  tension  have  been  ignored. 

From  a  mathematical  point  of  view,  the  identification  of  (13)  and  (16)  represents  only 
a  "first-order  match",  and  from  an  oceanographical  point  of  view,  it  leaves  us  without  a 
direct  "feel"  for  the  scales  of  the  motion.  Both  these  shortcomings  are  alleviated  by 
numerical  results^  for  the  beach  slope  c  «  ii/30.  Even  for  such  a  steep  beach,  there  is 
close  quantitative  agreement^  between  the  first-order  approximation  common  to  both  theories 
and  the  exact,  standing-wave  solution  of  the  classical,  linear  theory  over  a  distance 
extending  from  about  one-half  to  one  deep-sea  wavelength  from  shore.  The  length  of  this 
overlap  Interval  corresponds  to  about  one  local  wavelength,  at  such  a  large  beach  slope, 
and  at  the  deep  end  of  the  interval,  the  water  depth  already  exceeds  2/3  of  the  wave- 
penetration  depth.  The  values  of  u  a  in  the  overlap  interval  are  of  order  1/e,  and  the 
first  asymptotic  approximation  (Section  IV)  is  then  very  close  to  the  exact,  standing-wave 
solution  (Section  IV)  of  the  beach  theory.  In  short,  even  at  larger  beach  slopes  than  are 
common  in  nature,  there  is  a  substantial  interval  over  which  the  exact  solutions  of  the  two 
theories  differ  by  swich  less  than  could  be  detected  by  observation.  Remarkably,  moreover, 
the  beach  theory  is  seen  thus  to  describe  the  whole  shoaxlng  process  from  the  shoreline  to 
a  distance  normally  considered  as  one  at  which  the  water  is  "quite  deep";  in  particular,  it 
covers  all  distances  from  shore  over  which  the  water  would  be  considered  shallow  enough  for 
application  of  the  linear  longwave  theory. 

For  model-preserving  solutions,  the  "confidence  domain"  of  the  beach  equations  (1), 

(2)  is  therefore  very  much  larger  than  it  can  be  for  all  solutions  [F.  II].  It  extends  to 
distances  “  lOe  ^1^  from  shore,  at  least,  and  the  incidence  condition  [F.  (22)J  retains 
its  direct,  physical  meaning  over  (scaled)  times  t  of  order  e  \  at  least;  its 
mathematical  extension  [F.  IX]  is  dictated  only  by  the  requirement  of  time-uniformity 
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[F.  VIII]  and  by  the  objective  of  constructing  a  theory  not  at  all  limited  to  near-standing 
waves • 

A  different  aspect  of  these  issues  Is  Illuminated  by  the  question:  over  what  distance 
from  shore  do  even  "unbroken"  waves  display  s  seriously  "nonlinear"  character?  To  fix  the 

Vj 

Ideas,  if  2n  '  >  3.5  be  regarded  as  the  value  below  which  the  argument  In  the  Bessel 
functions  of  Section  IV  "is  not  large",  then  since 

(UgO/2)^  >  4u)^h*/(e^g)  -  8iilx*|/(eX_)  , 

asymptotic  approximation  of  the  Carrler-Greenspan  standing  wave  falls  only  for 

Ixgl  <  n,  |x*|  <  eX./2  . 

That  the  nonlinear  beach  theory  Is  needed  for  the  local  description  of  waves  over  only  such 
tiny  distances  from  shore  must  not,  of  course,  be  permitted  to  obscure  Its  qualitative 
long-range  influence:  it  Is  the  theorem  of  Section  III  wh.jh  makes  the  selection  between 
the  two  fundamental  solutions  of  the  linear  theory. 

The  critical  length  scale  L  is  defined  [F.  IV]  as  the  scale  of  distances  from  shore 
over  which  the  nonlinear  wave  mechanism  dominates.  As  a  scale,  of  course,  it  cannot  be 

assigned  a  unique  numerical  value,  but  any  reasonable  estimate  of  it  will  help  greatly  to 

Vo 

fix  the  ideas.  If  tti„o/2  »  2it  be  again  chosen  to  correspond  to  the  distance  from  shore 
B 

at  which  the  seaward  asymptotics  of  the  beach  theory  becomes  applicable  and  hence,  the  non¬ 
linear  aspect  of  that  theory  has  ceased  to  dominate,  then  the  distance  x  at  which 

u)  a/2  =  2  appears  a  reasonable  choice  of  the  critical  scale  L  for  illustrative  purposes. 
B 

It  corresponds  to  the  choice  L  =  i.e.,  o)^  ”  1.  Its  dimensional  value  is 

L  =  eX^/(2ii)  . 

To  illustrate  the  striking  implications  of  this  long-delayed  insight,  the  Table  lists 
these  values  of  L  for  various  beach  slopes  c  and  reference  wavelengths  X^. 


i'-  L’v.'.T'- '.’J'.'J  i  • 


■9 


spsunsn 


Table 


<*» 

1 

o 

1 

10“2 

10"^ 

to 

cn 

L  >  f.6  mm 

1 

m 

1.6  cm 

1.6  M 

0.16  imi 

10 

m 

16  cm 

1.6  cm 

1.6  mm 

100 

m 

16  cm 

1. 6  cm 

The  scale  L  Is  much  smaller  than  Intuition  could  have  suggested,  and  this  emphasizes 
again  the  dual  role  of  the  shore  singularity:  the  singular  enhancement  of  the  shoaling 
process  becomes  pronounced  only  in  the  immediate  neighborhood  of  the  moving  shoreline,  but 
the  beach  solutions  describing  it  also  describe  the  wave  motion  over  the  beach  at  much 
larger  distances  and  have  a  decisive  influence  on  the  wave  structure  in  the  deeper  sea. 

Ihe  small  values  of  L  in  the  table  suggest  that  surface  tension  should  play  an 
important  role  near  the  shoreline,  but  surface  curvature  there  depends  on  the  local 
amplitude  even  more  than  on  the  length  scale,  and  the  issue  Is  beyond  the  scope  of  the 
present  study. 
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VI.  Second  Friedrichs  Connection 

To  extend  such  considerations  of  a  match  between  shore  theory  and  linear  theory  on  a 
firm  basis  to  more  general,  developing  waves,  requires  a  useful  representation  of  such 
solutions  of  the  beach  equations  (1),  (2).  A  very  general,  Rlemann  representation  Is 
described  In  the  proofs  (Appendices  A,  B  of  [M] )  of  the  Incidence  and  Regularity  Theorems, 
but  that  representation  Is  complicated  due,  not  least,  to  Its  very  generality.  Waves  of 
oceanographical  Interest  Incident  from  the  sea  would  be  expected  to  be  quite  smooth,  and 
the  Incidence  Corollaries  tM.  4]  show  that  such  smoothness  carries  over  to  the  apparent 
solutions.  Such  waves  can  therefore  be  expected  to  possess  a  Fourier  transform  with 
respect  to  characteristic  time  X,  If  they  decay  sufficiently  as  X 

To  derive  the  Fourier  representation  frcm  Incidence  data  Is  complicated  and  can  be 
avoided  by  recalling  the  uniqueness  of  apparent  solutions  (M.  3].  If  Is  a  small, 

positive  number,  the  apparent  problem  Is  regular  for  a  >  o^,  and  In  Carrier's  notation 
(Section  IV),  the  solution  of  the  wave  equation  for  his  function  ^(0,X)  is  determined 
uniquely  for  o  >  0^  by  the  Cauchy  data  {♦(0^,X),  0^u(0^,X))  together  with  the  Initial 
data.  If  the  shore  conditions  IF.  (18),  (19))  are  to  loe  satisfied  as  well,  however,  the 
functions  ^(0^,X)  and  0^u(0|,X)  must  depend  on  each  other,  and  thus  ^(0j,X)  already 
determines  any  admissible  solution  uniquely. 

The  class  of  apparent  solutions  under  consideration  now  can  therefore  be  characterized 
by  the  assumption  that  ^(0^,X)  has  a  Fourier  transform 

8«g(«»)  “1^-0 

(since  ♦{0^,X)  i  0  for  0^  +  X  <  0  for  the  solutions  starting  from  rest)  in  the  strict 
sense,  if  ^(0^,X)  decays  sufficiently  as  X  *  »,  emd  In  the  distributional  sense,  in  any 
case.  The  non-dimensional  frequency  Is  here  again  ui  >>  Im^AI  In  terms  of  the  dimensional 
frequency  .  The  uniqueness  of  admissible  solutions  then  implies  that  It  must  have  the 


vVV 


tS-Cv.N 


Carrier  function 


because  this  solves  the  apparent  equation  and  has  the  correct  Fourier  transform  of 


'  V.-.-.  .V  ■  - 


This  confirms  the  plausible  expectation  that  any  beach  solution  of  amplitude 
small  enough  to  guarantee  model-preservation  can  be  represented  as  a  superposition  of 
standing  waves • 


It  now  follotfs  from  the  Initial  definition  of  ^  (Section  IV)  that 
1  2  1  3*  1  /*  .  iuo  -luX/2 

”8*7"  ■  7  ax  ■  TJI  /- 


U  -I4I.  ^  /-  ^  J,(^)e- 
o  3o  so  •'-»  B 


which  decays  even  faster  as  o  ♦  •  at  fixed  X,  so  that  the  Fourier  transform  of  n 


becomes  ultimately 


~  -i(D«g(Ul)  Jjj  (410/2) 

~  -2i(u/xo)^2J  (u,)  coa{iii0/2  -  x/4)  as  (u)o|  *  »  , 

0 

2 

and  from  [F.  (12)1,  also  X/2  ^  t  +  const.  Since  n_  decays,  moreover,  o  /16  -  -x  , 

B  B  B 

and  on  recalling  the  scaling  (11)  of  the  beach  theory,  010/2  ~  2oi„|ljc  j'^/u,  and  the 

Fourier  tremsform  of  n  is  approximated  by 
B 

~  (SOI  U)"  ^^2  |lx*|~  '''*co8{2ui,|Lx*  I'^^/U  -  s/4)  (18) 

B  B 

as  |oio|  *  o,  always  provided,  of  course,  that  the  Incident-wave  amplitude  Is  small 

enou^ . 


To  see  what  this  proviso  means  In  the  present  context,  note  that 


^  X  f°*  2.  ,  ,,  ,oia_ 

aX-  U  “  6g(4.)a^(— )do. 


and  therefore  from  (F.  (12)], 


at  1  _  au  _  -1  r, 

IX  2  2X1  J— 


201^6 

B  ,  ,410.,  -i4iX/2  doi 


Since  the  righthand  side  here  Is  1/2  for  5_(ai)  =  0,  it  has  a  positive  lower  bound  at 

B 

the  singular  line  o  »  0,  if  41^5^(41)  has  a  sufficiently  small,  upper  bound. 


The  more  general  wave  motions  now  under  consideration  possess  no  obvious  amplitude 


constant  and  accordingly,  also  no  obvious  definition  of  a  reflection  coefficient.  However, 
for  waves  of  shore  amplitude  spectrum  small  enough  to  preclude  "breaking",  the  inviscid 
beach  model  (1),  (2)  conserves  energy  and  therefore  implies  reflected  waves  of  amplitude 


comparable  to  that  of  the  Incident  waves.  The  waves  further  away  from  shore  can 


accordingly  be  expected  to  be  governed  by  the  linear  theory  to  a  satisfactory 
approximation.  To  make  that  theory  mathematically  determinate  requires  recourse  also  to 
Initial  conditions,  which  are  not  expected,  however,  to  be  very  relevant  to  what  one  wishes 
to  observe  in  practice.  It  will  therefore  be  more  helpful  to  assume  here  that  the  wave 
motion  far  from  shore  possesses  a  frequency  scale  fl  characteristic  of  its  spectrum,  which 
can  be  used  to  make  the  governing  equations  non-dimensional.  A  non-dimensional  time  t^, 
frequency  distance  Xj^,  surface  elevation  and  velocity  potential  may  then 

be  introduced  by 

•  *  2 

t^  =  Ot  ,  (iij^  =  iu,/0,  Xj^  =  -X  n  /g  , 

•  n*n^/q,  . 


In  view  of  the  expected  smoothness  of  the  waves  of  small  amplitude,  it  is  plausible  to 


assume  that  4^  also  possesses  a  Fourier  transform 


lidt 


in  the  same  conditional  sense  appealed  to  for  the  solution  of  the  beach  equations*  For 


every  fixed  (u  ,  ♦  must  then  satisfy  the  equations  solved  by  Friedrichs’*  and  hence* 

Li 


in  his  notation^,  where  Zp  =  Xp  +  iyp  is  referred  to  the  deep-sea  wave  number 
corresponding  to  i.e.. 


2  2*  2 
P  =  z  /g  =  OI^Z^ 

3 


The  Incident-wave  components  of  the  spectrum  are  then 


^  2CS_exp(-i(Zp-(i)t)  J 


in  the  deep-sea,  with 


C  =  -(21I€)'‘=  exp[-i»  /(8e)  +  iii/41 


so  that 

5„(0Ot  )  =  2C  S  (u.  ) 

represents  the  deep-sea  amplitude  spectrum  of  the  incident  waves,  and  of  the  reflected 
waves,  as  well,  because^  X.  *  0  in  the  deep  sea. 


and  In  that  approximation 

~  -  <_e“^'’)  •  (19) 

♦ 

When  a  compariaon  of  the  transforms  of  the  actual  surface  elevations  n  ^ 

«  2 

n  •  9^./^  predicted  by  the  two  theories  la  now  envisaged#  It  must  be  recognised  that, 
z* 

while  the  frequency  ratio  (ii./di  ••  V/lhSi)  la  fixed,  the  frequencies  themselves  range  over 

is 

all  the  real  numbers,  and  different  spectral  components  will  coaw  to  be  compared  at 
different  shore-distances  x*.  This  Is  not  meant  in  the  sense  that  the  real  beach  has 
uniform  slope  extending  to  Infinite  %«ater  depth  and  that  the  real  motion  remains  two- 
dimensional  over  such  Infinite  distances  from  shore*  Both  theories  admit  a  mathematical 
extrapolation  of  local  conditions  to  theoretical  conditions  at  other  locations  over  such  an 
Idealized  beach,  and  It  Is  this  extrapolation  which  provides  the  conceptual  clarification 
on  which  a  reasonably  simple  correlation  of  the  theories  can  be  based.  In  reference  to  a 
plane  beach  of  Indefinite  extent,  moreover,  a  comparison  of  different  spectral  components 
at  different  shore  distances  Is  natural  because  any  given  x  ^  0  corresponds  already  to 
effective  deep-sea  conditions  for  the  shortwave  components,  but  still  to  effective  shore- 
proximity,  for  the  longwave  components.  Given  any  ui^,  however.  It  was  seen  In  Section  V 
that  there  is  a  cooaon  range  of  x*  over  which  (18)  and  (19)  approximate  the  solutions  of 
the  respective  theories. 


-19- 


With  this  understanding,  a  comparison  of  the  approximation  for  given  by  (18) 

B 

•  2 

with  that  for  9h^/C2  given  by  (19)  becomes  meaningful  since  the  transforms  are  with 
respect  to  the  same  non-dimensional  time  in  the  limits  envisaged.  In  the  arguments  of  the 
trigonometric  terms,  |/U^  must  be  compared  with  u^x  /c,  and  they  are  identical 

Lt  la 

2 

because  =  co^/O  and  U  =  EgL,  by  definition.  The  trigonometric  terms  in  (18)  and 
(19)  are  therefore  identical,  if  and  only  if 

4^((Dj^)  s  -4_((iij^)  , 

by  which  the  comparison  is  seen  again  to  resolve  the  Indeterminacy  in  the  linear  theory. 

*  »  2 

The  identification  of  eLn_  and  gn,/fl  is  then  complete,  if  and  only  if 

B 

-  ieai,5  L^dtui.U)"  L~  s  -26_ga~^(iiu  , 


6(ia)  =  -2iii«  (di,  )U^/(eLn)^ 
o  “  L 

-  <’  +  lr>i 

2  -  L  Lfi^  4  2e 

v, 

with  u)  =  ui  n(L/eg)  This  determines  the  Fourier  transform  of  the  dimensional  solution 

JL 

of  the  beach  equations  uniquely  from  the  Fourier  transform  of  the  developing  wave  in  the 
deep  sea. 

The  relation  between  the  amplitude  spectra  5g(<»)  and  6^(u^)  given  by  the  last 
equation  still  depends  on  the  choice  of  the  scales  L  and  SI,  which  are  left  somewhat 
indefinite  by  the  mere  specifications  that  SI  is  a  scale  characteristic  of  the  deep-sea 
frequency  spectrum  and  L,  a  scale  characteristic  of  the  shore  distance  within  which 
"nonlinearity"  comes  to  dominate.  In  a  consistent  description,  the  two  scales  cannot  be 

completely  independent,  and  the  results  of  Section  V  show  that  L  must  be  envisaged  to  be 

2  2 
of  the  order  of  Eg/SJ  .  As  seen  at  the  end  of  Section  V,  moreover,  the  choice  L  =  Eg/O 

is  not  unreasonable  for  fixing  the  ideas;  it  ma)(es  m  *  and  simplifies  the  amplitude 


relation  to 


exp(^  (1  +  ^) 


Except#  perhaps,  at  the  longwave  end  of  the  wave  spectruin,  the  very  large  size  of  this 
ratio  is  an  unrealistic  feature  arising  from  the  neglect  of  dissipation  in  the  present, 


•  >  •  fc  ^  a  '  I 
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